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J^ . ABSTRACT 

The flux of cosmic-ray induced single hadrons near sea level has been measured 
with the large hadron calorimeter of the KASCADE experiment. The measurement 
corroborates former results obtained with detectors of smaller size if the enlarged veto 
of the 304 m^ calorimeter surface is encounted for. The program CORSIKA/QGS JET is 
used to compute the cosmic-ray flux above the atmosphere. Between Eq = 300 GeV and 
1 PeV the primary proton spectrum can be described with a power law parametrized as 
dJ/dEo = (0.15±0.03) ■Eq'^-'^^^^-^^ m''^ s''^ sr'^ TeY''^. In the TeV region the proton 
flux compares well with the results from recent measurements of direct experiments. 

Subject headings: cosmic rays - primary proton flux - air shower measurements 

1. Introduction 

The "knee" in the cosmic-ray energy spec- 
trum has been observed by many research 



groups and in several observables of air shower 
experiments. Typically, its position is found 
around a primary particle energy of 4 PeV. 
Many theoretical approaches to explain the 
knee exist. The most probable cause seems 
to be a superposition of spectra of many nu- 
clei each with an individual flux cut-off at dif- 
ferent energies. The hypotheses of the ori- 
gin and experimental findings, however, dif- 
fer significantly, in particular for the primary 
proton spectrum, which is of special rele- 
vance to understand cosmic-ray acceleration 
and propagation in the galaxy. A proton 
knee is claimed to be seen at different en- 
ergies, e.g. at 10 TeV by the MUBEE col- 
laboration (Zatsepin et al 1993), at 100 TeV 
by the Tibet group (Amenomori et al 2000) 
and at 4 PeV by the KASCADE collabo- 
ration (Ulrich et al 2001). When taking 
the experimental data of older direct mea- 
surements by balloon or satellite experiments 
above the atmosphere at face value, one might 
imagine a change of the power law slope at 
10 TeV. But recent measurements seem not 
to confirm such former conjectures. Precise 
measurements during the last years in the 
100 GeV region yielded proton fluxes lower by 
about 30% compared to older measurements, 
e.g. the new data of the BESS (Sanuki et 
al 2000), CAPRICE (Mocchiutti et al 2001), 
IMAX (Menu et al 2000) and AMS (Alcaraz 
et al 2000) collaborations. On the other hand, 
in the 100 TeV region recent publications on 
direct measurements report on higher flux val- 
ues, e.g. JACEE (Asakimori et al 1998) and 
Runjob (Apanasenko et al 2001). These find- 
ings indicate that the proton flux does not 
seem to decrease as strongly as anticipated so 
far. 

Therefore, it is of interest to determine the 
proton flux over a wide range of primary en- 
ergy using one single method. Such a method 
is the detection of single hadrons at ground 
level. These unaccompanied hadrons turn 



out to be intimately connected to primary 
protons. The latter penetrate deeper into 
the atmosphere with their hadronic compo- 
nent than heavier primaries of the same en- 
ergy and are the most abundant producers of 
single hadrons. Single hadron spectra have 
been measured using different experimental 
techniques like emulsion chambers, magnetic 
spectrometers, or calorimeters. In the past, 
measurements have been carried out both, at 
sea level (Cowan and Matthews 1971, Sio- 
han et al 1977, Fickle and Lamb 1979, Mielke 
et al 1994) and at mountain altitudes (Agli- 
etta et al 2003, Inoue et al 1997). Different 
definitions of single hadrons are used in the 
literature. For the present investigations a 
large calorimeter is used at sea level and sin- 
gle hadrons are defined as follows: Only one 
hadron with an energy of at least 100 GeV is 
reconstructed in the detector. In addition, the 
zenith angles are restricted to less than 30° in 
the analysis. During the last years the sim- 
ulation of air showers has improved consider- 
ably and the primary particle spectra can be 
deduced from ground based experiments with 
more confidence. 

The KASCADE calorimeter has been op- 
erating continuously and steadily for many 
years. Large data sets have been accumulated 
which allow to estimate the flux up to the PeV 
range. However, as air shower simulations in- 
dicate, more and more single hadrons turn out 
to originate from higher-mass primaries with 
increasing energy. Hence, the connections to 
primary protons become less stringent. 

2. Experimental apparatus 

The hadrons have been detected with the 
central calorimeter of the KASCADE exper- 
iment measuring cosmic rays near sea level 
(Engler et al 1999). It is a sampling calorime- 
ter consisting of layers of lead, iron, and con- 
crete absorbers interspersed with 9 layers of 
warm-liquid ionization chambers with an ac- 




Fig. 1. — Pattern of a single hadron event in 
the calorimeter. Each pixel represents the de- 
posited energy on an electrode of 25x25 cm^. 
The straight line indicates the reconstructed 
shower axis. 

ceptance area of 304 m^. The finely seg- 
mented chambers allow to measure the energy 
of individual hadrons and to reconstruct their 
point and angle of incidence. An example of 
a single hadron event is presented in figure 
1. The energy depositions in the ionization 
chambers are plotted from which a total en- 
ergy of 4.5 TeV for the hadron shown has been 
reconstructed. Apart from the hadron cas- 
cade no significant energy deposition is seen 
in the calorimeter layers. In particular, also 
the uppermost layer, where the entire electro- 
magnetic energy is measured, is nearly empty. 
The detection and reconstruction efficiencies 
of hadronic cascades in the calorimeter have 
been determined by simulating cascades with 
the detector simulation code GEANT (CERN 
1993)/FLUKA (Aarnio et al 1987, 1990). At 
100 GeV a trigger and reconstruction effi- 
ciency of 50% is reached, increasing to more 
than 70% at 500 GeV. Radiating muons can 
imitate a hadron. However, above 100 GeV 
the contributions of hadrons faked by muons 



is below the 1% level, for a detailed discus- 
sion see Mielke et al (1994). The maximum 
energy of an unaccompanied hadron detected 
was 50 TeV. During five years of operation 
(1996-2001) more than 1.5-10^ events have 
been recorded with at least one reconstructed 
hadron, out of which 2.5-10^ had one hadron 
only. 

The calorimeter is surrounded by an ar- 
ray of stations equipped with scintillators in 
which the electromagnetic and muonic com- 
ponents of an air shower are detected. A de- 
scription of the experiment can be found in 
Antoni et al (2003). 

3. Simulations 

The measurements have been accompanied 
by extensive shower simulations in order to 
understand the phenomenon of single hadrons 
and to determine the relation between the pri- 
mary proton spectrum and the single hadron 
spectrum at ground level. The program COR- 
SIKA 6.014 (Heck et al 1998) has been em- 
ployed with the code QGSJET 01 (Kalmykov 
and Ostapchenko 1993, Kalmykov et al 1997, 
Heck et al 2001) for high-energy hadronic in- 
teractions and GHEISHA (Fesefeldt 1985) for 
energies below 80 GeV. 

For the primary fluxes of nuclei above 
the atmosphere parametrizations according 
to power laws are taken from compilations 
by Wiebel-Sooth et al (1998) and Horandel 
(2003). In the latter review more recent mea- 
surements are taken into account and the 
parametrizations for proton, helium and iron 
have been updated. In total, about 2-10^'^ 
events have been simulated in the energy 
range from 100 GeV to 3 PeV for proton, he- 
lium, oxygen, and iron induced showers. The 
number corresponds to a data taking period 
of approximately 80 days for the calorime- 
ter acceptance. For illustration, the fluxes 
used are depicted in figure 2, marked with 
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Fig. 2. — Primary flux values for the indicated 
nuclei vs. the primary particle energy, accord- 
ing to parametrizations for oxygen (Wiebel- 
Sooth et al 1998) and for proton, helium, 
and iron (Horandel 2003). The correspond- 
ing spectral indices 7 are given. 

the corresponding spectral index 7 and ex- 
trapolated into the PeV region. The shower 
cores of the simulated events have been dis- 
tributed evenly over the calorimeter area ex- 
tended by 2 m at all four sides. The hadrons 
are tightly concentrated near the shower axis. 
The simulated distance distribution can be 
parametrized by an exponential function that 
falls off to 1/e within 4.3 m, nearly indepen- 
dent of energy. Using this lateral distribu- 
tion, calculations show that out of all primary 
particles, which are reconstructed as single 
hadrons 75% have been simulated. The simu- 
lations reveal, though, that for a given inter- 
val of hadron energy, the fraction of missing 
events is nearly independent of primary en- 
ergy. Since the missing percentage does not 
depend strongly on energy it has no signif- 
icant effect on the deduced primary proton 
spectrum which will be presented in section 5 
below. 
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Fig. 3. — Mean number of hadronic inter- 
actions in single hadron events vs. the pri- 
mary energy for three primary particles as in- 
dicated. 

Single hadrons belong to a particular class 
of air showers for which the detected hadrons 
stem from air showers with only a few inter- 
actions in the atmosphere. In vertical direc- 
tion the thickness of the atmosphere corre- 
sponds to about 11 interaction lengths. How 
many hadronic interactions occurred for sin- 
gle hadron events is shown in flgure 3. In 
the TeV range primary protons or helium nu- 
clei encounter 3 or 4 interactions only. In- 
tegrated over all relevant energies, the aver- 
age number of interactions is 3.6 for primary 
protons. Heavier nuclei do not generate sin- 
gle hadrons in this energy range, as revealed 
by the primary iron simulations. Oxygen and 
iron nuclei can be seen in this class of events 
only if the primaries have energies higher than 
30 TeV. 

The number of reconstructed single hadron 
events for the four classes of primaries are 
shown in figure 4 on the left-hand panel. On 
the right-hand the corresponding probability 
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Fig. 4. — Number of reconstructed single hadron events in the simulations vs. primary energy 
for the indicated nuclei (left-hand). The probability to find a single hadron vs. primary energy is 
shown in the right-hand graph. 



to find a single hadron event is plotted for 
primaries with respect to their energy. One 
observes that up to a few TeV, indeed, single 
hadrons originate mostly from primary pro- 
tons and that above 10 TeV also helium pri- 
maries contribute. At approximately 1 PeV, 
proton and helium contribute with equal num- 
bers to single hadrons and at higher energies 
also heavier nuclei have to be considered. The 
right-hand panel reveals that in the hundred 
TeV range 10% of all showers are of the sin- 
gle hadron type. This figure is valid for the 
present definition of a single hadron. If also 
accompanying electromagnetic energy would 
exclude single hadrons, they would be encoun- 
tered less frequently, as outlined in the next 
section. 

The simulations show how closely the single 
hadron events are related to the primary par- 
ticles. In figure 5 the distributions of primary 
proton energies are presented for three energy 
bins of single hadrons: for about 100 GeV, 



1 TeV and 10 TeV. The spread of the pri- 
mary energy is wide, but the bulk of parents 
have an energy ten times larger than the ob- 
served energy for all three intervals. Based on 
simulations it has been verified that the 25% 
of events missing do not change the shape of 
the distributions appreciably. Therefore, de- 
spite the large fluctuations, measuring the sin- 
gle hadron spectrum allows to deduce the pri- 
mary proton flux. 

4. Single hadron spectrum 

The flux of single hadrons is obtained using 
the trigger and reconstruction efficiencies as 
determined with the GEANT/FLUKA code. 
The data are given in table 1 and are pre- 
sented in figure 6. The errors quoted are es- 
timated systematic uncertainties concerning 
the fiducial area of the calorimeter, the to- 
tal data taking time, the effective solid an- 
gle, and the energy assignment, but are dom- 
inated by the trigger and reconstruction ef- 
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Fig. 5. — Frequency distributions of primary energies for three intervals of single hadron energies. 
Simulations with CORSIKA/QGSJET for proton initiated showers. 
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Fig. 6. — The single hadron spectrum: flux multiplied by the energy squared vs. single hadron 
energy. For comparison results from other experiments near sea level are presented. 
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Table 1: Single hadron fluxes from the vertical direction measured at sea level. The errors represent 
systematic uncertainties. 



ficiencies. They amount to approximately 
15% below 1 TeV, 25% in the TeV range and 
35% above 10 TeV. The numbers of collected 
hadrons given in the table indicate that sta- 
tistical errors can be neglected below 10 TeV. 
The data exhibit almost a power law in the 
double logarithmic graph. However, on the 
energy scale of over 2.5 orders of magnitude 
a gentle bend is apparent. Such a behaviour 
is observed in the simulations as well. For 
reference, measurements from the literature 
(Cowan and Matthews 1971, Fickle and Lamb 
1979) are shown as well as the data from the 
KASCADE prototype calorimeter (Mielke et 
al 1994). The last experiment exhibits some- 
what larger fluxes compared with the present 
data due to its smaller surface of 6 m^. The 
KASCADE calorimeter with 304 m^ fiducial 
area has a more efficient veto for multiple 
hadron detection. Especially for low hadron 
energies this may cause the differences. The 
data of the two older experiments, which both 
had smaller apertures of about 0.65 m^- sr, 
show a similar shape. 



In the present investigations also accom- 
panying electromagnetic radiation detected in 
the scintillators of the array stations can be 
accounted for. The graph in the left-hand 
panel of figure 7 shows the number of sta- 
tions which have registered at least one min- 
imum ionizing particle in coincidence with a 
single hadron. In fact, only 6% of the single 
hadron events are not accompanied by a sig- 
nal in the array stations. As can be inferred 
from the insert, a small probability exists that 
all 252 array stations have a signal in coinci- 
dence with a single hadron. All this signi- 
fies that the notion "single hadron" is some- 
how artificial, it depends on the experimen- 
tal conditions and ipso facto changes from ex- 
periment to experiment. Demanding that not 
more than eight stations have an electronic 
signal, i.e. ^dep > 5 MeV, already reduces 
the number of single hadron events by 1/3. 
This can be seen in the right-hand panel of 
figure 7, where the single hadron spectrum 
is compared with this additional requirement. 
One observes that the electron veto becomes 
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Fig. 7. — Left-hand panel: Number of electromagnetic detectors with ^dep ^ 5 MeV. Right-hand 
panel: The single hadron spectrum compared with additional vetos: Not more than eight array 
stations or no station have registered a minimum ionizing particle. The ordinate has been multiplied 
by the energy squared. 



stronger for large hadron energies which in 
the mean originate from higher primary en- 
ergies. Simulations indicate that the sensitiv- 
ity to primary protons is also enhanced. In 
the figure the flux with no signal at all in the 
252 stations is presented as well. However, 
because the number of events fulfilling this 
requirement is too small, the veto condition 
of not more than eight stations activated has 
been chosen for further analyses, in particular 
when deducing the primary proton spectrum. 

5. Primary proton spectrum 

The single hadron spectrum is converted to 
a flux of primary protons by attributing to the 
measured hadrons a probable primary energy 
according to the energy distributions which 
are illustrated in figure 5. In principle the 
single hadron spectrum 5'(lgEHad) has to be 
converted to a flux of primary protons /(Ig Eg) 
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Fig. 8. — Spectral index of the primary proton 
spectrum deduced from the measurements vs. 
the slope used in the simulations of the prob- 
ability density distributions. 
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Table 2: Primary Hux of protons deduced from the single hadron spectrum assuming protons only 
(second column) and with contributions from helium and heavy nuclei subtracted (third column). 
The errors are systematic uncertainties (in parentheses the contribution from the systematic errors 
of the single hadron spectrum). 



by solving the integral equation 



^(IgEHad) = / A(lgEHadllgEo) /(IgEo) dlgEg 

(1) 
A(lg Euadl Ig Eq) is the transfer function trans- 
forming the primary flux spectrum into the 
measured single hadron spectrum at ground 
level accounting for the interactions in the 
atmosphere. Several methods exist to de- 
convolute one-dimensional spectra. For sake 
of simplicity a slightly different approach has 
been chosen, which turned out to be rather ro- 
bust and straightforward. Knowing the prob- 
ability density distribution i?(lgEo| IgEnad) 
for a given single hadron energy IgEjjad) the 
primary proton flux can be inferred by 



/(IgEo) = y 5(lgEo|lgEHad)5(lgEHad)dlgEHad. 

(2) 
Using the calculated probability distribu- 
tions (see e.g. figure 5) assuming primary 
fluxes as given in figure 2, the primary proton 
flux is obtained. The resulting energy spec- 
trum depends only slightly on the slope of the 
a priori flux assumptions, as can be seen in 
figure 8. Plotted are the power law indices de- 
rived from the data versus the assumed indices 



in the simulations of the probability distribu- 
tions. One realizes that the method yields 
rather stable results. 

The outcome for the uncorrected primary 
fluxes are given in the second column of ta- 
ble 2, assuming that all primaries are pro- 
tons. They are multiplied by E^'^ in order 
to compensate the steeply falling power law 
spectrum. As seen in figure 4, at higher ener- 
gies also helium and, to some extent, the heav- 
ier nuclei initiate single hadron events. These 
contributions are subtracted by taking into 
account the parametrized primary fluxes as 
shown in figure 2. The final proton fluxes are 
presented in the third column. In the TeV re- 
gion the difference of the flux values amounts 
to about 20%, a value comparable to the er- 
rors as discussed in the following. 

The errors quoted are estimates of system- 
atic uncertainties. They include the calcu- 
lation of the probability density distribution 
and the errors in the single hadron spectrum 
(given in parentheses), both added quadrat- 
ically. The uncertainties in i?(lgEo| IgEnad) 
are estimated to amount to 15% below a few 
TeV and 20% at higher values of Eq. 

Graphic representations of the results are 
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Fig. 9. — The flux of primary protons as function of energy (black and grey shaded squares). The 
dashed line represents a fit of a power law. The lines indicate the maximum systematic errors 
estimated. For comparison results of direct measurements are shown: AMS (Alcaraz et al 2000), 
BESS (Sanuki et al 2000), CAPRICE 94 (Boezio et al 1999), CAPRICE 98 (Mocchiutti et al 2001), 
HEAT (DuVernois et al 2001), Kawamura (et al 1989), IMAX (Menn et al 2000), JACEE (Asakimori 
et al 1998), LEAP (Seo et al 1991), MASS (Bellotti et al 1999), MUBEE (Zatsepin et al 1993), 
RICH (Diehl et al 2003), Runjob (Apanasenko et al 2001), Ryan (et al 1972), Smith (et al 1973), 
Sokol (Ivanenko et al 1993), and Webber (et al 1987). 



shown in figures 9 and 10. The present 
data are plotted as filled squares and the 
maximal errors are indicated by the two 
lines. The error bars on the individual 
points represent the systematic uncertain- 
ties in the single hadron fiux. The squares 
follow a power law reasonably well, a cor- 
responding fit yields dJ/dE^ = (0.15 ± 
0.03) • ^-2-78±o.03 ^-2 g-i g^-i XeV-i which 

is indicated by the dashed line. It should be 
kept in mind that for higher energies above 
100 TeV contributions of helium and heavy 
nuclei of up to 50% had to be subtracted. 
In the figures the corresponding values are 
marked as shaded points. 

In figure 9 the data are compared to results 
of direct measurements above the atmosphere 



taken from the literature. One recognizes dif- 
ferences between the individual results of the 
order of a factor of two. Within these uncer- 
tainties the present data are compatible with 
results from the literature. At energies around 
100 GeV the most recent data of direct mea- 
surements scatter at the lower bound of the 
published fluxes. Our proton fiux extrapo- 
lates well to these data. This fact can be inter- 
preted in such a way that in this energy region 
from 10^ to 10"^ GeV the hadronic shower cas- 
cade within the atmosphere is well described 
by the program CORSIKA with the interac- 
tion code QGSJETOl. Also other tests have 
proven that below 1 PeV the latter code de- 
scribes the shower propagation best (Antoni 
et al 1999 and 2001). 
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Fig. 10. — Flux of primary protons as presented in figure 9. For comparison results of indirect 
measurements by air shower are shown: EAS-TOP (Aglietta et al 2003), Tibet (Amenomori et 
al 2000), HEGRA (Aharonian et al 1999), and SYS (Inoue et al 1997). 



In figure 10 our data are shown together 
with fluxes for primary protons obtained by 
experiments using indirect methods of mea- 
surements as well. Within the errors given, 
the data corroborate previous measurements. 

6. Summary 

Using the large hadron calorimeter of the 
KASCADE experiment during three years of 
effective data taking, 2.5T0^ events have been 
accumulated for which a single hadron was 
reconstructed. With these data the energy 
spectrum of single hadrons has been derived. 
These data are somewhat lower than fluxes 
published previously. This is attributed to 
the large surface area of the calorimeter which 
acts as a more efficient veto against multi- 
hadron events. 

Single hadron events are particular air 
showers, which predominantly stem from pri- 
mary protons in the energy region considered. 



Applying large sets of simulated single hadron 
events, and assuming primary particle fluxes 
as obtained by direct measurements extrapo- 
lated into the PeV region by power laws, the 
primary spectrum for protons could be de- 
termined. The data extend over nearly four 
orders of magnitude in primary energy and 
can be described by a power law with a flux 
parametrized as dJ/dEo = (0.15 ± 0.03) • 
£;-2.78±o.03 ^-2 g-i gj.-i XeV-^ In the TeV 

region the proton fluxes agree well with recent 
measurements of direct experiments above the 
atmosphere. 
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